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We consider a version of the extended Clifford Group which is defined in terms 
of a finite Galois field in odd prime power dimension. We show that Neuhauser's 
result, that with the appropriate choice of phases the standard (or metaplec- 
tic) representation of the discrete symplectic group is faithful (as opposed to 
merely projective), also holds for the anti-unitary operators of the extended 
group. We also improve on Neuhauser's result by giving explicit formulae 
which enable one easily to calculate the ( ant i-) unitary corresponding to an ar- 
bitrary (anti-)symplectic matrix. We then go on to find the eigenvalues of an 
arbitrary (anti-)symplectic matrix. The fact that in prime power dimension 
the matrix elements belong to a field means that this can be done using the 
same techniques which are used to find the eigenvalues of a matrix defined over 
the reals — including the use of an extension field (the analogue of the complex 
numbers) when the eigenvalues are not in the base field. With the eigenval- 
ues of the (anti-) symplectic matrix in hand it is straightforward to find the 
eigenvalues, the order and all roots of the corresponding (anti-)unitary. We 
then give an application of these results to SIC-POVMs (symmetric informa- 
tionally complete positive operator valued measures). We show that in prime 
dimension our results can be used to find a natural basis for the eigenspace 
of the Zauner unitary in which SIC-fiducials are expected to lie. We go on 
to use the extension field to construct a parameterization of the displacement 
operators analagous to the one commonly employed in quantum optics, using 
the complex eigenvalues of the annihilation operator. Finally, we apply our 
results to the MUB cycling problem. Wootters and Sussman have shown that 
in every even prime power dimension there is a single Clifford unitary which 
cycles through a full set of MUBs (mutually unbiassed bases) . We show that 
in odd prime power dimension d, although there is no Clifford unitary, there 
is a Clifford anti-unitary which cycles through the full set of Wootters-Fields 
MUBs if d = 3 (mod 4). Also, irrespective of whether c? = 1 or 3 (mod 4), 
the Wootters-Fields MUBs split into two groups of {d + l)/2 bases in such a 
way that there is a single Clifford unitary which cycles through each group 
separately. 
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1. Introduction 



Since its original introduction into the field of quantum information the Clifford 

group [iiiiiaiiiaiBiHiiiiiniiiiiiiiiiiiiiiiaiMiiiiiaiiaisniiss^ 

[24l [25l [26l [27l [28l [29l [30l has found numerous applications. Our own interest in the 
group, and the extended Chfford Group ^12,^ of which it is a subgroup, is due to 
their role in the theory of mutually unbiassed bases (MUBs) pl l ITS l IT9 l [2 H [25 l [26 l 
EF, '3T, 'M", 'M', 'M', '3F, 'W, 'FT, "M", 'W, '411 ia 1^ 

[55t ,54. ,55. ,56, ,57, 58,, ,59,, ,60. 61, 62., 63^ and symmetric, informationally complete 
positive operator valued measures (SIC-POVMs) [ia[ni[2lllIlll3[5Ql[511[55l[5S 
leOlIMIMlIMliZllMlIMlIZSIZllI^IZaiiaiZ^ However, the 

results we derive along the way may be of wider interest. 

We should remark at the outset that it is possible to identify at least three 
different (though closely related) mathematical constructs which commonly go by 
the name "Clifford group" . We will clarify the precise sense in which wc arc using 
the term in Section [2] (see paragraph following Eq. ([28])). For now wc will confine 
ourselves to saying that we are concerned with a version of the Clifford group which 
is defined over a finite field 32j, and which therefore only exists in prime power 
dimension. There are some significant complications when the dimension is a power 
of 2, so in this paper we confine ourselves to the case of odd prime power dimension. 
We hope to treat the even prime power case in a subsequent paper. 

In the first part of the paper (sections [2] to [5]) we prove a number of general 
results concerning the Clifford and extended Clifford groups. In the second part 
we apply these results to a problem which arises in the theory of SIC-POVMs and 
to the MUB cycling problem. We also show how they can be used to construct 
an alternative parameterization of the displacement operators. In a subsequent 
paper we will use them to complete the proof that the minimum uncertainty states 
introduced by Wootters and Sussman [25l |2^ and Appleby, Dang and Fuchs [77] 
exist in every prime power dimension. 

The Clifford group is built out of two kinds of operator: displacement operators 
Du, labelled by vectors of the form 



with components in the finite field {d being the Hilbert space dimension), and 
symplectic unitaries Up, labelled by 2 x 2 matrices of the form 



with elements in and determinant = 1. The set of all such matrices constitutes 
the (discrete) symplectic group, SL(2,Frf). The symplectic unitaries permute the 
displacement operators according to the prescription 



Explicit expressions for the operators and Up will be given below. 

The extended Clifford group is built out of the above operators together with 
anti-symplectic anti- unitaries Up, for which the matrix F has determinant = — 1. 
We refer to the set of all 2 x 2 matrices with either sign of the determinant as the 
extended symplectic group ESL(2,Fd). 




(1) 




(2) 



UpD^Uj, = D 



Fu 



(3) 



2 



After reviewing what is already known on this subject in Section [2] we begin 
our analysis in Section [3] by showing that the phases of the operators Up can be 
chosen in such a way that the map F ^ Up (the metaplectic representation [HIH]) 
becomes a faithful representation of the group ESL(2,Fd), as opposed to one that 
is merely projective. In other words 



for all Fi, i^2- The fact that this is possible for the group SL(2,F(i) was shown 
by Neuhauser 4]. We establish that the result extends to ESL(2,Fd). We also 
improve on Neuhauser's result by giving explicit formulae which enable one easily 
to calculate Up for any given matrix F. In fact we give two different sets of formulae 
for, in addition to giving explicit expressions for the matrix elements, we show (in 
Section!!]) that U p can be simply expressed as a linear combination of displacement 
operators. 

We next (Section [5]) consider the problem of finding the order of an arbitrary 
symplectic/anti-symplectic operator Up {i.e. the smallest positive integer m such 
that [/™ = 1), its eigenvalues and eigenspaces, and its roots. To solve this problem 
we make essential use of the fact that is a field. This means that we can use 
exactly the same techniques to diagonalize a matrix 



with elements € that we would if its elements were real numbers. In particular 
we can use the trick of going to an extension field. If it should happen that the 
characteristic equation of a matrix over the reals has no real solutions, then we 
simply embed the reals in the field of complex numbers. Similarly here: if the 
matrix F has no eigenvalues in F^; we can embed F^ in the extension field Vj2 . Once 
one knows the eigenvalues, finding the order is a comparatively simple matter. We 
give explicit formulae which enable one easily to calculate the order of an arbitrary 
matrix F S ESL(2,Fd). With these formulae in hand it is straightforward to 
calculate the order, eigenspaces and eigenvalues, and roots of the corresponding 
unitary/anti-unitary Up. 

In Section [6] we apply the methods developed in the previous sections to a prob- 
lem which arises in theory of SIC-POVMs. Scott and Grassl's [50] recent exhaustive 
numerical investigation encourages the speculation that, in all dimensions, every 
Weyl-Heisenberg covariant SIC-fiducial is an eigenvector of a canonical order 3 
Clifford unitary [12] . Generally speaking this unitary belongs to a different version 
of the Clifford group than the one considered here. However, if the dimension is 
prime (as opposed to a power > 1 of a prime) the two versions are the same, and 
the results obtained in this paper become applicable. We show how they can be 
used to find a natural basis for the eigenspace of the order 3 symplectic unitary in 
which SIC-fiducials are expected to lie. 

In Section [7] we describe an alternative labelling of the displacement operators. 
In quantum optics it is usual to parameterise the (continuous variable) displacement 
operators with the complex variable 



Up^Up^ — Up^p^ 



(4) 




(5) 




(6) 
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(where q, p are the quadratures). A similar construction can be carried through in 
the discrete case, using the extension field F(i2 instead of C. Taking the discrete 
logarithm we obtain an integer labelling of the displacement operators, which leads 
to a natural way of representing ESL(2,Frf) as a group of permutation matrices. 
This result is needed in Section [TOl but it may be of some independent interest. 

Finally, in Sections [8] to [10] we consider the MUB-cycling problem. Wootters 
and Sussman [251 I26j have shown that in every dimension equal to a power of 
2 there exists a single Clifford unitary which cycles through a full set of MUBs. 
We investigate the situation in odd prime power dimension. Wc show that in this 
case there is no single Clifford unitary which cycles through the full set of Wootters- 
Fields MUBs [33j (at least for the version of the Clifford group considered here — see 
below). However, if d = 3 (mod 4) there is a single anti-CliSoid anti-umtary which 
cycles through them. Furthermore, it is possible, for all odd prime power d, to split 
the Wootters-Fields MUBs into two groups of {d + l)/2 bases each in such a way 
that there is a single Clifford unitary which cycles through each group separately 
(we say that a unitary with this property is half-cycling). This leads to a natural 
labelling scheme for the MUBs, in which a cycling anti- unitary (when it exists) 
increases the integer index by 1 and a half-cycling unitary increases it by 2. 

Wootters and Sussman [25] used the existence of cycling unitaries in even prime 
power dimension to prove the existence of minimum uncertainty states in all such 
dimensions. Sussman [26] subsequently extended the proof to an infinite subset 
of the prime power dimensions = 3 (mod 4). In a subsequent paper we will use 
the results obtained in Sections [8] to [TOj of this paper to complete the proof, and to 
show that minimum uncertainty states exist in all prime power dimensions, without 
exception. 



We begin by reviewing some relevant definitions and known facts concerning 
Galoisian variants of the Clifford and extended Clifford groups [H [H [SI [SI [13 [HI 



[ni[Ta[20l|2l|22|23[2i[2i[3l[MllMlliaiMllZll^ This will also 



allow us to clarify the relation between the version of the Clifford group considered 
in this paper and definitions used by other authors. 

Let us begin by defining the Galoisian variant of the Weyl-Heisenberg group 
(or generalized Pauli group). The ordinary Weyl-Heisenberg group is defined by 
choosing an orthonormal basis |0), . . . , |d — 1) in d-dimensional Hilbert space 
and then defining operators X and Z by 



where u = e~ and addition of indices is mod d. If d is odd the Weyl-Heisenberg 
group is the group generated by these operators (in even dimension there is a slight 
complication — see, for example, ref. [12j). In dimension d — p"^, where p is a prime 
number and n is a positive integer, the Galoisian variant of the Weyl-Heisenberg 
group is defined similarly except that instead of labelling the orthonormal basis 
by the integers mod d one labels them by the elements of the finite field (for 
a summary account of the aspects of the theory of finite fields which are relevant 
here see, for example, Vourdas [20], Pittenger and Rubin [38] or Klimov et al [54] : 
for a more comprehensive treatment see, for example, Lidl and Niederreiter [82]). 




2. Preliminaries 



X\x) = + 
Z\x) = uj'-'\x) 



(7) 

(8) 
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One then defines, for all x, m e F^, 

Xu\x) = \x + u) (9) 
Z^\x) ^ u'<^^^\x) (10) 

where w = e~ (observe that uj is now a p*"^ root of unity, not a d^^ root of unity 
as in the case of the ordinary variant), and where tr : — > Zp is the field theoretic 
trace defined by 

n-l 

tr(x) = ^ x^P") (11) 

r=0 

(making the natural identification of Zj,, the integers modulo p, with the set {z € 
F^: = z}). Notice that if n > 1 we cannot write = X^, = because it 
does, in general, make no sense to raise a Hilbert space operator to the power of an 
element of a finite field. However if n = 1 we can identify F^ with Zp, which means 
we can make these replacements, and the Galoisian variant of the Weyl-Heisenberg 
group becomes identical with the ordinary variant. 

Next, for each vector u ~ (^1,1*2) G F^, define the displacement operator 

^u = t"-("i"^)X„,Z„, (12) 

where r — Note that t'^ — uj and that r, like u, is a p^^ root of unity. We 

have, as an immediate consequence of this definition, 

Di = (13) 

and 

D^D^ = t''^'^^ D^+^ (14) 
where (u, v) is the symplcctic form 

(u, v) — tr(u2Ui — U1V2) (15) 



It follows that 



(Du)'' = Dkr (16) 



for every integer k and all u G F^. In particular 

{D^f = 1 (17) 

So the displacement operators are all order p (apart from Dq). 

It can be seen from the above that the set of operators {w^Du : u e F^, m e Zp} 
constitutes a group, which is what we are calling the Galoisian variant of the Weyl- 
Heisenberg group, and which we will denote Wd- 

Many authors define in a way that may look, on the face of it, rather different, 
as a tensor product of n copies of the ordinary variant of the Weyl-Heisenberg 
group in dimension p (the "many-particle" definition, in Gross's [T71 [571 ISI] 
terminology). To see that this definition is in fact equivalent to ours let be any 
basis for the field F^, and let be the dual basis. So 

tr(eres) = 6r,s (18) 
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for all r, s. For arbitrary x G¥d let Xr ~ tr(a;er) (respectively Xr — tr(a;er)) be its 
expansion coefficients relative to the basis Cr (respectively Sr). So 

n n 

X — ^^^^ Xj^Gj' — ^^^^ XfGf (X9) 
r=l r=l 

Let Hp (respectively Hd) be p-dimensional (respectively rf-dimcnsional) Hilbert 
space, and let |0), . . . ,\p—l) be the standard basis for Tip. Let D^^ be the ordinary 
Weyl-Heisenberg displacement operators in dimension p. So 

i:>P|fc) =T''i"^+2'=''^|fc + Ui) (20) 

Finally, let S : Hd Hp (X" . . . Hp be the linear isomorphism whose action on the 
basis states is 

S\x) = \xi) ■ ■ ■ \x„) (21) 

Then 

= i?K,„,) = ® • • • ® ^(..„.a.„)) S (22) 

for all u (where uir = tr(?iie,.), U2r = tr(?i2e,.)). This defines an isomorphism of the 
groups Wp (E" • • • (8) Wp and Wd- Note, however, that it is not a natural isomorphism 
(because there are many pairs of dual bases for the field F^^). 

These definitions (the one we gave earlier, and the one which consists in iden- 
tifying Wd with the tensor product Wp (8) • • • (g) Wp) both have their advantages. 
The tensor product definition is, perhaps, the more natural of the two in a situa- 
tion such as occurs when one has a register consisting of n different p-dits. In a 
situation like that the tensor product structure reflects the actual physics. Also the 
tensor product definition may be found attractive because it makes no reference to 
the mathematics of Galois fields, which is unfamiliar to many physicists. However, 
those advantages partly disappear if one is interested in a single qudit which doesn't 
naturally split into a product of n different p-dits. In that case the imposition of an 
(arbitrary) tensor product structure introduces a needless complication. Moreover, 
even in a situation where the system of interest does in fact consist of n distinct 
p-dimensional subsystems one might want to make use of the properties of Galois 
fields. In a case like that, given that one is going to have to introduce them even- 
tually, it would seem to make sense to introduce Galois fields right at the outset, 
in the definition itself. In this paper we are going to make very heavy use of Galois 
fields. We consequently prefer our definition. 

Let us now turn to the definition of the Galoisian variant of the Clifi^ord group. 
Here again we have a choice between two altc;rnative definitions. However, whereas 
in the case of Wd it is essentially a matter of taste which definition one prefers, 
here the choice is substantive in that the groups defined are non-isomorphic. Let us 
begin with what we will call the unrestricted (Galoisian) Clifford group, Cd- This 
is the set of unitaries U which reshuffle and re-phase the displacement operators, 
according to the prescription 

UD^U^ = e'9(")Z?j(„) (23) 

for all u and suitable tZ-dependent functions f,g. Since D^^ = 1 for all u the phase 
e»5(p) must be a p^^ root of unity. So we must in fact have 

f/£>uf/^ =w''("^-D/(u) (24) 
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for some function h taking values in Zp. In other words Cd is the normahzer of Wd 
{i.e. it is the set of unitaries U such that UWdU^ ~ Wd). 
It follows from Eq. (fH|) that, for any such U, we must have 

for all u, V. So the function / must have the property 

/(u + v) =/(u) + /(v) (26) 

In particular 

/(mu) = /(u+... + u) = m/(u) (27) 

for all TO e Zp. We will refer to a function with this property as p-linear. 

We define the restricted (Galoisian) Clifford group to be the subgroup of Cd 
comprising those unitaries which satisfy the additional requirement that / be, not 
merely p-linear, but d-linear in the sense 

J{xn) = xj{n) (28) 

for aU X e Frf. 

Some authors (for example Cormick et al [21]) use the term "Clifford Group" to 
mean Cd] others (for example Vourdas [21]) use it to mean CJ. Gross [TT ] [22 l [27| . IST] 
has proposed to resolve the ambiguity by referring to Cd as the "many-particle" 
Clifford group, and to as the "single-particle" Clifford group. We have not 
adopted Gross's terminology ourselves because it seems to us that the mathematical 
distinction between the groups Cd and CJj has nothing specially to do with the 
distinction between composite and non-composite physical systems. Suppose we 
have a non-composite physical system, and suppose we adopt the "single-particle" 
definition (in Gross's terminology) of Wd, as represented by Eqs. (HJ) and (fTOjl . Then 
we will still be faced with the fact that the normalizer of Wd is C^, not C^. So Cd 
is no less relevant to single-particle systems than it is to many-particle ones. The 
reason is that the distinction between p-linearity and d-linearity does not depend 
on the assumption of a tensor product structure. It therefore seems to us that our 
terminology, referring to Cd as the unrestricted group and to as the restricted 
group, is more appropriate. 

In this paper we will focus on the restricted group, C^. This is not because 
we consider to be of greater physical and/or mathematical importance than 
Cd- Nor is it because we are more interested in single-particle systems than many- 
particle ones. It is just that the d-linear transformations of have a number of 
nice properties which make them worthy of separate investigation. In particular, 
they have all the properties needed to prove the existence of minimum uncertainty 
states in every odd prime power dimension (which was the problem that originally 
motivated this research). We have therefore chosen to focus on in this paper, 
and to defer a discussion of Cd to what we hope will be another, subsequent paper. 

If it is assumed that / is d-linear Eq. (|25)) implies 

(a) The function / is of the form 

/(u) = Fu (29) 
where F e SL(2,Fd) — i.e. is a 2 x 2 matrix 
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with entries in ¥d and such that det F = 1. 
(b) The function g is of the form 

g{u) = {x,n) (31) 

for some fixed vector x ^ i^d)'^- 
We refer to the matrices in SL(2,Fd) as symplectic matrices. 

Let us sketch the proof of these statements. It is easily seen that (-Fu, Fv) = 
tr(det F{u2Vi — ^1^2))- In view of Eq. (^5)1 this means 

2h{u + v) - 2h{u) - 2h{v) = tr((det F - 1)(m2Wi - U1V2)) (32) 
(mod p) . Interchanging u and v we also have 

2h{u + v) - 2h{u) - 2h{v) = ~tr((detF - l)(?i2Wi - U1W2)) (33) 

(mod p). It follows that detF = 1. Setting detF = 1 in Eq. (f32|) we deduce that 
h is p-linear. In terms of the dual bases introduced earlier we have 

n 

U = (wi,U2) = ^(Mlr.(er,0) +U2r(0, e^)) (34) 
r=l 

where Ujr — ti{ujer). So if we define 

(n n \ 

-^/i((0,e,))e,,^/i((e,,0))eJ (35) 
r=l r=l ) 

it follows from the p-linearity of h that 

/i(u) = (x,u) (36) 

This establishes that for every C/ e there is a matrix F g SL(2, F,^) and vector 
X e (Frf)2 such that 

JJDuC/^ ^c^^^'^^L'fu (37) 

for all u. The converse is also true: for each F and x there is a corresponding 
U e C\. In fact, it can be shown that for each F G SL{2,¥d) there is a unitary Up 
such that 

UpD^Ul = (38) 

(in Section [3] we will give an explicit formula for Up)- So the unitary U — UpF)^ 
will transform the displacement operators according to the prescription of Eq. ([37]). 
Moreover it is easily seen, in view of the irreducibility of the standard representation 
of Wd and Schur's lemma, that if U' is any other unitary satisfying Eq. ([57]) then 
U' = e^^U for some phase e*^. In short consists of the set of unitaries of the 
form 

e^'^UpD^ (39) 

for arbitrary F e SL(2,Fd), x G ^d)'^^ G M. The properties of the displacement 
operators are easily understood, so in the remainder of this paper most of the work 
will go into understanding the properties of the symplectic unitaries Up ■ 

It is also interesting, and for our purposes important, to consider anti-unitary 
operators U with the property 

UD^U^ = e''^"^-D/(u) (40) 
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for all u and suitable functions /, g. As before the fact that = 1 means that 
e«9(") must be a p**^ root of unity. So we must have 

C/i?uC/^ -w"(")i?/(u) (41) 

We will refer to such operators as anti-Clifford anti-unitaries. The group which 
consists of the Clifford unitaries £ Cd together with the anti-Clifford anti-unitaries 
just defined we will call the extended (Galoisian) Clifford group 8d- 

The description of £d parallels the discussion in Section 4 of ref. [12]. Let J be 
the matrix 

-0 

(so det J = — 1) and let Uj be the anti-unitary which acts by complex conjugation 
in the standard basis: 



(43) 



Then it is easily seen that 

UjDJj] = Dj^ (44) 
for all u. Now consider the unitary operator V = UUj. We have 

implying that V E Cd- Conversely, ii V E Cd then VUj is an anti-Clifford anti- 
unitary. So £d is the disjoint union 

Sd^CdU {CdUj) (46) 

In this paper we will be concerned, not with the full group £d, but with the 
subgroup 

£:^=C5u(Q[/j) (47) 

We accordingly define FiSL{2,¥d) to be the group consisting of all 2 x 2 matrices 
with entries in F^; and determinant = ±1. We refer to the matrices in ESL{2,¥d) 
with determinant = — 1 as anti-symplectic matrices. If detF = 1 let Up he the 
unitary defined earlier, while if detF = — 1 let it be the anti- unitary 

Uf - UfjUj (48) 

Then it is easily seen that £^ consists of all operators of the form 

U ^ e'-'^UpD^ (49) 

with F e ESL(2,Fd), X £ (Fd)^ and e*^ an arbitrary phase. If detF = 1 then U is 
a Clifford unitary; if deti<" = — 1 then it is an anti-Clifford anti- unitary. In either 
case 

UD^U^ ^lo'^^^'^'^Dfu (50) 

for all u. 
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3. Faithful Representation of ESL(2,F(i) 

Let Fi,F2 be any two elements of ESL(2, ¥d). Then 

UpUpD^Ulul^ = ^F.F.u = Up^pD^Ul^^^ (51) 

for all u. It follows from the irreducibility of the standard representation of Wd 
and Schur's lemma that 

UF.UF.^e^'UF.F, (52) 

for some phase e*^, meaning that the map F Up is a projective representation 
of the group ESL(2,F(i) (the metaplectic representation [31[52). In this section we 
show that it is possible to choose the phases of the operators Up in such a way that 
the representation becomes, not merely projective, but faithful, so that 

UfiUf^ = Uf^f^ (53) 

for all Fi, F2. The fact that this is possible for the unitaries corresponding to 
the elements of SL(2,Fd) has been shown by Neuhauser [4]. We improve on his 
result by (a) showing that it extends to the whole of the group ESL(2,Fd) and 
(b) giving explicit formulae which enable one easily to calculate the unitary/anti- 
unitary corresponding to an arbitrary element of ESL(2,Fd). 

We begin by establishing some preliminary results. Let F^ (respectively Z*) be 
the set of non-zero elements of F^ (respectively Zp). Let be a primitive element 
for the field F^, and for each a: S FJ let logg x be the unique integer in the interval 
[0, d—2] with the property 

X = e^°^o ^ (54) 

Let Q (respectively N) denote the set of quadratic residues (respectively non- 
quadratic residues) of F^: 

Q ~ {x <e¥*i^ : X = for some y G F^} (55) 
N = {x e¥*^:xT^y^ for aU y e F^} (56) 

Clearly x ^ Q (respectively x G N) if and only if logg x is even (respectively odd). 
Consequently Q and N each contain exactly {d — l)/2 elements. Let l{x) be the 
quadratic character of F^: 

fl xeQ 

l{x) = < -1 xeN (57) 

[0 X = 

and let lp{z) be the quadratic character of "Lp {i.e. the Legendre symbol (^)). Note 
that /, like lp, is multiplicative: l{xy) = l{x)l{y) for all x,y G F^. We then have the 
following result 

Lemma 1. Suppose z e Z*. Then 

\lp{z) 

Proof. Define 

e = = (59) 



if d is an even power of p , , 

(58) 

if d is an odd power of p 
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Then = e, implying that e G Zp. Moreover e'' = 1 if and only if r = (mod 
p — 1), implying that e has multiplicative order p — 1. So e is a primitive element 
for Tip. By construction 

logeZ= (l+p+---+p"-i)log,z (60) 

for all z e Z*. If n is even then {1 + p+ ■ ■ ■ is even, implying that logg z is 

even for all z € Z*. Consequently l{z) ~ 1 for all z G Z*. If, on the other hand, n 
is odd then {1 + p + ■ ■ ■ + p"^^) is odd, implying that logg z is even if and only if 
logj z is even. Consequently l(z) — lp(z) for all z e Z*. □ 

It is also convenient to define the quantity 

(61) 

II X = 

It has the property 
Lemma 2. For all x £ ¥d 

{lix)y^li-x) (62) 

Proof. Suppose d is an even power of p. Then it follows from Lemma [T] that 
Z(— 1) = 1. Consequently l{—x) = l{—l)l{x) = l{x). Moreover, the fact that n and 

?i(p + 3) 

p + 3 are both even means that i 2 is real. The claim is now immediate. 

Suppose, on the other hand, that d is an odd power of p. Then in view of 
Lemma [T] and a standard result of number theory (see, for example. Hardy and 
Wright [83]) we have 

l{-x) = lp{~l)l{x) = i-l)^l{x) = (-l)-^^/(a;) = i-"(P+3);(x) (63) 

The claim now follows. □ 

We will also need the following Gaussian sum formula 

Lemma 3. For all a e¥*^ and b E ¥d 

J2 r^r(a.-+»-) = ^/dl{a)T- (64) 

Remark. The fact that we work in terms of r rather than u> means that we can give 
a single formula which covers both the cases p = ±1 (mod 4). 

Proof. Using Theorem 5.15 in Lidl and Niederreiter [82] in conjunction with the 
fact that '}2,x&a ^^'^'^^'^ = we find 

= l{a) 1{x)lu''^''^ (66) 
j {-iy'-^l{a)Vd p = 1 (mod 4) 



(-l)"-ii"/(a)\/d P= -1 (mod 4) 



(67) 
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Since {d+ l)/2 is the multiplicative inverse of 2 considered as an element of Zp we 
have r*'''^"^ ' = w*''^"^ and, consequently, 

^tr(a.^) ^ U-ir-'l{a/2)Vd p = 1 (mod 4) 

ij'/ \{-ir-'i"l{a/2)Vd p=-l(mod4) 

Suppose, now, that n is even. Then it follows from Lemma [T] that l{a/2) — 1(a). 
Moreover 

, f(-l)"-^^(«) p=l(mod4) 
^"'^ [(-l)"-ii"/(a) p = -l(mod4) ^ 

Consequently 

irrespective of the value of -p. 

Suppose, on the other hand, that n is odd. Then it follows from Lemma [T] and 
a standard result of number theory (see, for example. Hardy and Wright 83 ) that 

Z(a/2) = ;,(2)Z(a) = |^("^ P = ±1 (mod 8) 

' \-l{a) p = ±3(mod8) ^ ^ 

So 

'l{a)Vd p = 1 (mod 8) 

^tr{ax') ^ i -i''l{a)Vd p^3 (mod 8) 



-l{a)Vd p = 5 (mod 

J"l{a)^/d p = 7 (mod 



Comparing with the definition of I (a) we again find 

r*^('^-') = Vdlia) (73) 

irrespective of the value of p. 
Finally, if 6 7^ 0, we have 



□ 



We are now ready to establish the main result of this section. Let 
be any matrix e ESL(2, Frf). If det F = 1 define 

If, on the other hand, det F — —I define 

Uf = UfjUj (77) 
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where Ufj is the unitary defined by Eq. (I76|) . J is the anti-symplectic defined by 
Eq. (|42l), and C/j is the anti-unitary defined by Eq. ^ (note that det(FJ) = 1, 
so this definition makes sense). 
We have 

Theorem 4. For all F E ESL(2,Fc;) the operator Up is unitary ifdetF = 1 and 
anti-unitary ifdetF = — 1. We have 

(1) For all F e ESL(2, F<j) and ueWj 

UpDJjl - (78) 

(2) For all F, F' G ESL(2,Fd) 

UfUf' = Uff' (79) 

= (80) 
Proof. Suppose, first of all, that detF = 1 and (3 — 0. Then 

UfD^uI^ J2 r''^"'''-^'-y"^\x\D^\y)\ax){ay\ (81) 

^ r*'-(""^(-'-^'))T"-("i"^+2""«^2')5,,^+„„Ja;)(2/| (82) 

^ ^tr(aui(7«i+5«2))^tr((7Ui+5«2)!y)|y (83) 

x,ye¥d 

= Df^ (84) 

which establishes Eq. ([78)l for this case. Setting u = we also see that Uf is 
unitary. 

Suppose, next, that detF = 1 and f3 y^O. Then 

Xi,X2, 



i ^tr(/3 i(a(yJ-y^)-2(a:iai-a:2a2)+<S(a:?-x^))+«i«2+2«2y2) 

d ^ 



Xi,X2, 
Vl,V2&fd 



^ ^yi,y2+ui {x2\ 



2:i,2:2eFd y ^ y2eFi 

^ ^tr(/3"^(i5(xi-a:2)-2ui2;i+Q«i+pjiiii2jj|2.^^^2;2| 
— tr(/3-i(5(a:?-x^)-2«i2:i+Q«?+/3Mi«2))x , ,^ InWrnl 

£Ci,a;2eF£j 

^ Y ^tr((ani+/3«2)(7«l+A«2))^tr((7«i+5«2)a;)|2, _^ _^ PU2) {x\ 

xeWd 

= Df^ (85) 

which establishes Eq. ((78)) and the fact that Uf '^s unitary for this case also. 
It remains to consider the case when det F = —1. We have 

UjD^ul = Dj^ (86) 
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In view of the foregoing it follows that 

UpD^ul = UpjUjDJj]ulj = (87) 
We now turn to the proof that the representation is faithful. Let 



F=\Z f) (89) 



be two matrices G ESL(2,Fd), and let 

be the product F1F2. 

To start with suppose that det Fi — det F2 — I. There are four cases to consider 

Case 1. /3i = /32 = 0. Then 

= l{a) J2 T'''^°'^'''^\ax){x\ 
xe¥d 

= Up (90) 
Case 2. (3i ^ Q, (32 ^ 0. Then 

Ez 

Vd 

l{-ail32) ^ ^tr(a-i/3-i((7i/32+5i52)2;=^-2a;y+aia2?/=)) l^^^^l 



Xi ,X2 , 

y2e¥d 



= Up (91) 
Case 3. /3i ^ 0, P2 = 0. The proof is similar to that of Case 2. 
Case 4. (3i ^ 0, /32 7^ O.Then 

^ K-Pl)K-P2) ^tx[l3-^{atyl-2xiyi+Stxl)+li-\a2yl-2x2V2+S2xl)) 

Xl,X2, 

vi,y2&¥d 

X 5y,,:,2\xi){y2\ 

^ l{~Pl)l{-P2) ^ I^^tr{l3-^p-^l3z''--2(l3-^x+l3-^y)z)\ 

X ^trif3^'S,x'+l32''^2y')\x)(^y\ (92) 

It follows from Lemma [3] that 

^ ^ti{0^^l32^P^^-'^i0i^x+f32^y)z) 
26Fd 
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So if /3 = 

Uf.Uf, = r(-/3i)r(-/32) M^^^^'<-^^''^+^^''^^^y') \ - f3,f3^'y){y\ (94) 
To evaluate this expression observe that if /3 = 



ai Pi\ _ fa \ / ^2 ~f32\ _ f aS2 -a/?2 



(95) 



from which it is easily seen that 

a = (96) 
7 = -P^'^i - Pi'c^2 (97) 

Also, in view of Lemma [21 

f(-/3i)/(-/32) = !{-l3i){l{f32)y = ?(-/?i/32) = /(-/3i/32~') = ^(«) (98) 
Substituting these expressions into Eq. we deduce 

C/;^,[/j^, = l{a) J2 r*'("^^')|mj>(y| = C/;^ (99) 
Suppose, on the other hand, that P — 0. Then Eq. (|92p becomes 

[/f,c/i., = -^/(-/?i)U-/?2)U/3rV2~'/?) 

X ^ /r(/3-i(/32-^(/3a2-/3i)y^-2a;y+ft-H/35i-/32)2^'))|2.^^y| (IQO) 

To evaluate this expression observe that 

f3a2 -Pi= 132a (101) 
P5i - /?2 = /3i<5 (102) 

Also 

i{-PSK-P2)KPi^p2'P) = ^(-/3i)Z(-/32)(^'(-/3r'/32"V))* 

= -l-^l{-Pl)l{-P2)l{-Pl'P2'P) 

?i(p + 3) 

^-z-^;(-/3) 

= /(-/?) (103) 
Substituting these expressions into Eq. (|100p we deduce 

We now turn to the case when either or both of -Fi, F2 is an anti-symplectic 
matrix. It is readily verified that for all F g ESL(2,Fd) 

UfUj = Ufj (105) 

UjUf = UjF (106) 
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Now let Fi be symplectic and F2 anti-symplectic. Then in view of the foregoing 

Uf.Uf, = Uf,Uf,jUj = Uf,f,jUj = Uf,f, (107) 
It then follows from this and Eq. (|106p that if Fi is anti-symplectic 

UfiUf2 — Ufi,jUjUf2 ~ UfijUjf-2 = UF1F2 (108) 
for all F2, symplectic or anti-symplectic. Finally, we note that Eq. ([80| is an 
immediate consequence of Eq. ()79|) and the fact that the Uf are all unitary/ anti- 
unitary. □ 

4. Expressing the Uf in Terms of the £)„ 

In the last section we gave explicit expressions for the matrix elements of the 
operators Uf in the standard basis. However, it is sometimes useful to express 
them instead as linear combinations of the The fact that this is possible is an 
immediate consequence of the fact that 

Tr (dId^) = d<5u,v (109) 

which shows that, relative to the Hilbert-Schmidt inner-product {A, B) — Tr(A^^_B), 
the operators ■^^'u are an orthonormal basis for the space of linear operators on d- 
dimensional Hilbert space (note that we use the symbol Tr for the ordinary matrix 
trace, and tr for the field- theoretic trace). The following theorem gives the relevant 
formulae. 



Theorem 5. Let 

-J 

he any matrix € SL(2, F^) (so det F = 1 ). Let t = Tr(F). Then 
(a) 



(110) 



whe 



(b) Ift 



F (112) 



2-t 



Ere.. r~ "(^'^^^^0,,. (3 = and J ^0 (113) 

.Do /3 = 7 = 

We also have 

(c) Ift ^2 

(d) Ift = 2 

( l{-p)Vd p^O 

Tt{Uf) = I l{j)Vd p = andj^O (115) 

Id (3 = j = 
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Proof. Suppose that /? = 0. Then it follows from the definitions of I?u, Up that 

TT(DiUp) = {^\Di\y){y\UF\x) 



l{a) r*'-("^^'-2"-^+"^"^)4,,-.i'5,,a. 

\/d[(7)T-t'-('^"'"2)5„^ i = 2and77^0 (116) 

(i(5u,o t ^2 and 7 = 

where we used Lemma |3] in the last line. Suppose, on the other hand, that /3 ^ 0. 
Then 

Tr(z)tc/^) = {ADi\y){y\Up\x) 

— ^tr(^P^'-{ax^-2xy+Sy^)-2u2y+uiU2) g 

x,yeVd 

^ ^tr{l3-\t-2)x^+2{p-\S-l)ui-U2)x+t3-^5ul-uiU2) 

^(l{t-2)T(--^-) t^2 
\\/d/'(-/3)T*'^(''""')<5„,^/3-i(i_,)„, t = 2 

where we again used Lemma |3l Statements (c) and (d) of the theorem are special 
cases of these formulae (since Do = 1). Statements (a) and (b) are immediate 
consequences of them and the fact 

□ 

5. Eigenvalues and Orders of the (anti-)Symplectic Matrices 

We now turn to the problem of calculating the eigenvalues and order of an 
arbitrary symplectic/anti-symplectic matrix 

Once these are known it is straightforward to find the order, the eigenvalues and 
eigenspaces and the roots of the corresponding unitary/anti-unitary Up- 

Let t = a + 5 he the trace (in the ordinary matrix sense) of F, and let A be its 
determinant. We say that F is 

(1) type 1 if - 4A e Q 

(2) type 2 if ^2 _ 4A e TV 

(3) type 3 if t2 _ 4A = 0. 

Suppose, to begin with, that F is type 1. Then F is conjugate to the matrix [23l 

EH El [85] 
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This matrix is diagonalizable. In fact 

dct(F - XI) = -t\ + A (121) 
So the equation det(F — XI) — has the two solutions 



A. = (122) 



where y/f^ — 4 A is one of the two elements of with the property (\/i^ — 4A) = 

— 4A (guaranteed to exist because of the assumption that t"^ — 4A e Q) and 
where, as usual, division by 2 means multiplication by the multiplicative inverse of 
2 considered as an element of F^. It is easily verified that 

A = A+A_ (123) 

t = X+ + X- (124) 

This means we can write 

A+ = r (125) 

A_ = AO-'' (126) 

where r = log^A+. The fact that (A6'-'' - r)^ = t'^ - AA G Q means that 
AO-"- 7^ 9\ So the matrix 

(127) 



Ar 1 

is a well-defined element of ESL(2,F£i). It is straightforward to confirm that 

So the problem of finding the order of F reduces to the problem of finding the 
smallest positive integer m such that 



AO-'') ~ [q 1 



(129) 



If A = 1 we immediately deduce that 



where [r, (d — 1)] is the greatest common divisor of r and d — 1. Suppose, on the 
other hand, that A = — 1. Then the order of F is the smallest even positive integer 
m such that 0™ = 1. So 

if J^JJ^ is even 
o^djF) = { (rf-i) . (131) 

We can write this formula more compactly if we observe 



d-1 



'l^'d- 1] if is even 

Ur^d-l] if is odd 



(132) 
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from which it follows 

Let us note that for all r in the interval < r < d — 1 the matrix 

^-(? Or^to-r) (134) 

is a well-defined element of ESL(2,Frf). However, for some values of r the matrix 
is type 3 and not diagonalizable. It is easily seen that the values of r for which F 
is a type 1 matrix with eigenvalues O'' , A9~^ are 

1 < r < d- 1, r 7^ ^^-^1 ifA = l (135) 

l<r<d-l, ~ 1) I if A = -landd=l (mod4) (136) 

{r: 1 <r < d-1} if A = -1 and d = 3 (mod 4) (137) 

In particular the maximum order of a type 1 matrix is d — 1, the maximum being 
achieved by (for example) matrices for which t — 6 + A6~^. 

Suppose, next, that F is type 2. Then F is again conjugate to the matrix [231 

[281 Ell ES] 

-f) (13a, 

However the fact that the discriminant — 4A G N means that the equation 

det(F - A/) - a + A = (139) 

has no solutions in the field F^. We can deal with this problem in the same way that 
we would if we had a matrix defined over the reals whose characteristic equation 
had a negative discriminant: namely, we can go to a suitable extension field, in this 
case F(i2 . Let ^ be a primitive element for ¥^2 . We make the identification 

¥d^{xe Fd2 : x'^-'^ = 1} (140) 

F|J thus consists of all powers of It follows that d"^^^ is a primitive element 

for ¥d which, without loss of generality, we can identify with 9. 

To define the square root of an arbitrary element x € N observe that x — ^'^(''+1) 
for some odd integer k. So if we define yG; = 9 2 we will have x = (y/x) . 
Observe that, since 9 ^ = — 1 and k is odd, 

{^y^9^^9^^-^ (141) 

Now let 

A± ^— (142) 

be the solutions of Eq. (|139p . Using the fact that (x±yY = x'^±y'^ for aU x,y e ¥^,2 
we deduce 



4 = ^ (143) 



In view of Eq. (|141|) and the fact that x'^ = x for all a; G F^ it follows 

Xl - A^ (144) 
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So 

A+ = e'' (145) 
A_ = e*"^ (146) 

for some integer k. Since A+A_ = A, and taking account of the fact that 9 ~ = 
— 1, we must have 

mod {(P — 1). So if we define rj = 

A+ = (148) 

A_ = r]'-'^ = (-1) V (149) 

for some integer r which is in the range 1 < r < 2(d+l), and which is even if A = 1 
and odd if A = —1. 

Since — 4A ^ the eigenvalues are distinct and so the matrix 



is a well-defined element of ESL(2,F(j2). It is readily verified that 

5F5-i=f^ . (151) 



(-l)V 

from which we see that the order of F is the smallest positive integer m such that 
j^mr _ Since 77*^ = 1 if and only if = (mod 2{d + 1)) we conclude 

It follows from this that every symplectic matrix has order < d+ 1, and every 
anti-symplectic matrix has order < 2d + 2. Let us note that there exist both 
symplectic and anti-symplectic matrices which achieve these upper bounds. This 
is a consequence of the more general fact that for every value of r in the interval 
1 < r < 2{d + 1) there is a matrix F e ESL(2,Fd) with eigenvalues rj'', (-l)''ry~''. 
To prove this observe 

{f]'^ + i-iyr^-'^y = 7]*' + {-IYt^-'''^ = + {-IYti-'' (153) 

(since r]'^^ = {—lYri~^), implying that r]^ + {—lyrj^^ G F^. So the matrix 

^=f? .ri'/Cl-rl (154) 



^1 ,,'- + (-1)-,?-'- 

is in ESL(2,Fd) and has eigenvalues rf, (— l)'"r?~''. In particular 



-1 T + 'n 

is a symplectic matrix having the maximum order d+1 and 

^0 1 

1 v~ 'n~^ 

is an anti-symplectic matrix having the maximum order 2(rf -|- 1). 



(156) 
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It remains to consider type 3 matrices. A type 3 matrix is not diagonalizable 
except in the trivial case when F — ±/. We have the further complication that the 
conjugacy class of a type 3 matrix 

F=(: (157) 



v7 ^ / 

depends, not only on the values of t and A, but also on the values of f3 and 7. 
Finally, we need to consider the cases d — I (mod 4) and d = 3 (mod 4) separately. 
Using the results in ref. [25] we give in Table [T] a representative for each of the 12 
conjugacy classes when d = 1 (mod 4), while in Table[2]we give a representative for 
each of the 4 conjugacy classes when d — 3 (mod 4). The last column in each table 
gives the order of the matrices in the corresponding conjugacy class, calculated 
using the formula 



X y \ fx'' rx^ 



Ox \ X 



(158) 



We have thus proved 



Theorem 6. Let F be any element of F,SL{2,¥d). Let t be the trace of F (in the 
ordinary matrix sense) and let A be the determinant. Then 
(1) IfF is type 1 let 




-t + y/f^ - 4A \ , , 

r^logJ ^ j (159) 

Then the order of F is given by 

''-^ «/A = l 

o'-d(^^) - r z/A^-1 ^^^^^ 

where [u,v] denotes the greatest common divisor of u and v. The allowed 
values of r are 

r: 1 < r < d- 1, r 7^ ^^-^1 «/ A = 1 (161) 

r: 1 <r < d-1, r^ ^^'^^ if A ^ -1 and d ^ 1 (mod i) (162) 

{r:l<r<d-l} tf A ^ -I and d = 3 (mod 4) (163) 

In particular the maximum order of a type 1 matrix is d — 1, the maximum 
being achieved by (for example) symplectic matrices for which t = 9 + 9^^ 
and anti- symplectic matrices for which t — — 9~^ . 
(2) If F is type 2 let 9 be a primitive element ofV^p and let t] = 9^~ . Define 



-t + Vt^ - 4A \ , , 

r^logj ^ j (164) 

Then the order of F is given by 

Every value of r in the interval 1 < r < 2{d + 1) is possible, r is even for 
symplectic and odd for anti- symplectic matrices. So the maximum order of 
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a type 2 symplectic matrix is d + I, the maximum being achieved by (for 
example) matrices for which t = rf' + rf''^ , while the maximum order for a 
type 2 anti- symplectic matrix is 2{d + 1), the maximum being achieved by 
(for example ) matrices for which t = rj + rj'^ . 
(3) If F is type 3 the order is as tabulated in Tablel^ (if d — 1 (mod A)) or 
Table\M (if d ^ 3 (mod 4)). 



A 


t 


A 7 


representative 


order 


1 


2 


/3 = 7 = 


(I 0\ 
\ 1 / 


1 


1 


2 


/3 or 7 e Q 


V 1 


P 


1 


2 


13 or J £ N 


(X .\ 
I 1 / 


P 


1 


-2 


= 7 - 


/-I 0\ 


2 


1 


-2 


/? or 7 e g 


( -\) 


2p 


1 


-2 


P or ^ e N 


( -0 


2p 


-1 


2i 


/3 = 7 = 


A \ 


4 


-1 


2i 


/3 or 7 e g 


i\ 


4p 


-1 


2i 


(3 or J e N 


/ i J/ \ 

* / 


4p 


-1 


-2i 


/3 = 7 = 


(-i \ 
lo -^) 


4 


-1 


-2i 


/3 or 7 e g 


(o ) 


4p 


-1 


-2i 


/? or 7 e 




4p 



Table 1. Orders of Type 3 matrices when d — \ (mod 4). Here 
V is any fixed element of A^, and i is one of the two solutions to 
the equation = — 1 (so i does not have its usual meaning — for 
instance if d = 5 then i = 2 or 4). Note, also, that it is not possible 
for /? € g and 7 e iV, or for /3 G iV and 7 G g. For more details 
see ref. [28j. 



With these results in hand it is straightforward to find the order, eigenvalues and 
eigenspaces and roots of the unitary/anti- unitary Up corresponding to an arbitrary 
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A 


t 


/?, 7 


representative 


order 


1 


2 


= 7 - 




1 


1 


2 


/3 or 7 7^ 


(i 


P 


1 


-2 


= 7 = 


( -"0 


2 


1 


-2 


/3 or 7 7^ 


( -0 


2p 



Table 2. Orders of Type 3 matrices when d — 3 (mod 4). Note 
that for these values of d there are no type 3 anti-symplectics. For 
more details see ref. [28j. 



matrix F G ESL(2,F(i). In fact, let m = ord(F). Then m is also the order of Up- 
So the eigenvalues of Up are , with r = 0, 1, . . . , m — 1. The projector onto the 
eigenspace of Up with eigenvalue is 

Tn—l 

Pr = — Ve"^t/^ (166) 

s=0 

The dimension of the eigenspace (possibly zero) is Tr(Pr)- Finally, suppose we want 
to find the s'*^ roots of Up- i-e. the set of (anti)-symplectic unitaries Uc (possibly 
empty) such that Uq = Up. Clearly Ua has this property if and only if = F. 
Suppose, for the sake of definiteness, that F is type 1, so that 

for some S, r. Then G is an s*'' root if and only if 

for some t, A such that st = r { mod d — 1) and A" = detP. The problem of 
finding the roots of a type 2 or 3 matrix is equally straightforward. 

6. SIC-POVMs 

The results in the last section are potentially relevant to the problem of con- 
structing SIC-POVMs (symmetric informationally complete positive operator val- 
ued measures) in prime dimension. Every known [5D] SIC-POVM which is covariant 
under the Weyl-Heisenberg group can be constructed from a fiducial vector which 
is an eigenvector of a symplectic unitary Up for which Tr(_F) = — 1 (mod d). It 
is therefore of some interest to find a natural basis for the appropriate eigenspace. 
The results in the last section enable us to do that in every odd prime dimension 
(since if d is prime the Galoisian variant of the Weyl-Heisenberg Group considered 
in this paper coincides with the ordinary variant) . 
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Throughout this section it wiU be assumed without comment that d is prime 
number > 5 (our analysis does not apply to the case d — 3) . li t = Tr(F) = — 1 
and F is symplectic then — 4A = —3. It is a standard result that —3 G Q 
(respectively —3 G N) if d = 1 (mod 6) (respectively c? = 5 (mod 6)) (see, for 
example, Hardy and Wright [83]). So F is type 1 if d = 1 (mod 6) and type 2 if 
d = 5 (mod 6). 

Consider first the case d = 6m +1 = 1 (mod 6). Without loss of generality we 
can take F to be the matrix 



l2m 







Define 



- ( ^?-2™ ) (169) 



Then Up ~ C/q™. Uq is the permutation matrix 

d-l 

UG = -J2\er){r\ (171) 

Define 

^ S — 1 

where a = and r = 0, 1, . . . , 6m — 1. Also define 

IV^) = |0) (173) 
Then it is readily confirmed 

Ug\A)='j1A) (174) 

for r = 0, 1, . . . , d — 2, and 

Uam = m (175) 

Consequently 

Uf\A) = AHVr) (176) 
for r = 0, 1 , . . . , 6m — 1 (where A — e ^ ) , and 

= 1^^) (177) 

It appears from Scott and Grassl's exhaustive numerical investigation |80j that the 
fiducial lies in the eigenspace of highest dimension, which is the one with eigenvalue 
1 spanned by the 2m + 1 vectors |V'o), \tpo), IV'a), • ■ • , IV'em-a)- 

Now consider the case d = 6m —1 = 5 (mod 6) . Without loss of generality we 
can take F to be the matrix 

F=(^5 :}) (178) 

(the Zauner matrix [65]). This is a type 2 matrix, and therefore not diagonalizable. 
However, we can still use the techniques described in the last section to find an order 
6m symplectic matrix G such that F = G^™. It will turn out that the eigenvalues 
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of the unitary Ug are all non-degenerate, which means that it provides us with a 
natural basis for the subspace in which the fiducial is hypothesized to lie. 
To see this observe 




^=(4^ 'J^-']{ 'n .-4™) I ',;-^] (179) 



So if we define 
where 




(180) 



2r _ „ 2r 

7r = -i ^- 181 

the matrix G will have the desired property F — G^™ . Note that we can use the 
same trick to calculate an arbitrary power of G: 




^-4„ 
4r7i_ 



(182) 



Note also that, for all r. 



( 2r -2r\£; 2dr _ -2dr 

d ^ [V ) ^ ^ V ^ (183) 

^jj2m rq—2m'^d j^2dni j^—2dm ^ ^ 

(since rj'^ = —ri~^). So even though it is defined in terms of the quantity 77, which 
lies outside the field F^, 7r itself lies in the field F^. Finally, it should be observed 
that the matrix G is not unique. In fact if k is any integer relatively prime to 
(P — 1 then r is also a primitive element of ¥^2 . So 77 in the above formulae can 
be replaced with Ty'^, for any such k. For the convenience of the reader one possible 
choice of the matrix G is listed in Table [3] for each of the first 18 primes = 5 (mod 
6). 

Now consider the unitary Ug- Since Ug order d+1 its eigenvalues are powers 
of cr = . The projector onto the eigenspace of Ug with eigenvalue is 

6m — 1 

= 7^ E ^""'^G (184) 

s=Q 

To find the dimensions of the eigenspaces we need to calculate Tr{UQ) = Ti-{Ug=)- 
This can be done using Theorem [5l Let tg = Tr(G''). We have 

ts-2 = 7„+, + j„i-s = iv" ~ y^'f (185) 

So = 2 if and only if s = (mod Qm). Also 

{7f - f]-'f = ry"^" - ry-'^" = (-l)'+'(??' - V'l (186) 



25 



d 


G 


d 


G 


d 


G 


5 






47 1 
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) 
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y 4z 
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00 Q ^ 

— ZZ o 
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107 1 


,-31 -16, 


) 


17 


-0 




59 


/ 9^ 9R\ 

1,-28 3 J 




113 


1,47 9 ) 




23 


rn -7^ 

^7 4y 




71 ( 


^ 8 34 ' 
^-34 -29, 


) 


131 1 


^-44 11 ' 
,-11 -33, 


) 


29 


fi2 -9^ 

^9 3y 




83 


/18 -25\ 
1,25 -7) 




137 1 


^-58 20 ' 
,-20 -38, 


) 


« ( 


' 7 If 

-11 18 


) 


89 1 


^-25 -12' 
,12 -37, 


) 


149 


f47 -12\ 
1,12 35 ) 





Table 3. One possible choice for the matrix G for each of the first 
18 prime dimensions = 5 (mod 6). Note that other choices are 
possible. 



(since rj'^" — {—lyij So if s 7^ (mod 6m) then y/tg — 2 G if and only if s is 
odd. Consequently 

l{t, - 2) = (-1)"+^ - -cr^"" (187) 
for all s 7^ (mod 6m). We also have 7s = if and only if s = (mod 3m). So it 
follows from Theorem [5] that 

Tr(C/^) = -cr^™^ (188) 
if s 7^ (mod 3to) . It is easily seen that 

G^™=("o^ (189) 
Taking into account the standard result [83] 

Z(-l) = (-l)^=-(-l)™ (190) 
we deduce that Eq. (|188p continues to hold when s — 3m. Consequently 

1 / '^""^ \ 

^'^^'■^ " 6^ r " ^ = 1 - Sr,3m (191) 

So the eigenspaces are all one dimensional, with the exception of the eigenspace 
corresponding to r = 3m which is zero dimensional. 

Let {tpr) be the normalized eigenvector corresponding to the eigenvalue cr'" (notice 
that \ipr) can be calculated, up to a phase, from the explicit formula we have given 
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for Pr). We have UF\tpr) = X^ltpr), where A — ct^™ — . Let Sk be the eigenspace 
of Up with eigenvalue A'^. Then bases and dimensions of these subspaces are as in 
Table |31 Scott and Grassl's exhaustive numerical investigation [5(T suggests that 
fiducial vectors always exist in iSi and iS2, and that if to = (mod 3) they also exist 
in So- 



subspace 


basis 


dimension 


So 


IV'o), IV'3>, • ■ • , IV'3m-3), IV'Sm+a), • • ■ , IV'em-s) 


2to- 1 


Si 


\lp4), |V'6m-2) 


2m 


S2 


\i'2), \lp5), • ■ • , |V'6m-l> 


2m 



Table 4. 



7. Representing the F Matrices as Permutation Matrices 

The matrices F e ESL(2,F) permute the vectors u e (F^)^ which label the 
displacement operators. For many purposes it would be convenient if we could label 
the vectors by an integer r, so that F acts according to the rule F: — > Ufi^^) 
for some function / acting on the index (it will, for instance, be convenient in 
Section [TU|) . Of course one could assign the integer label arbitrarily. However, one 
would like to find a natural way of doing it, such that the function / takes a simple 
form for every F e ESL(2,Fd). 

We can achieve this by borrowing an idea from quantum optics. In the analysis 
of continuous variable systems it is often convenient to replace the two quadratures 
q, p with the single complex variable 

a = ^{q + ip) (192) 

One can obtain a discrete analogue of this by letting the extension field ¥^2 play 
the role of C. We accordingly define, for each u e (F^)^, 

x{u) = rj^^ui + U2 (193) 

(one is free to use any basis for ¥^2 in this definition; we choose the basis {?7~^, 1} 
because that will prove convenient in Section [TU]) . The map u x{u.) is easily seen 
to be a bijection of (F^)^ onto ¥^12, with inverse 




X ^ u{x) = \ I (194) 

This suggests that we define 

Dl = (195) 
( "e" for "extension field labelling" ) . We then have 

Z?S^^'=T<-'"'>i?S+.' (196) 
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for all X, x' e , where 



X'^^x' — XX' 



(^^^0=tr r;^^:t (197) 



Let us now consider the action of an arbitrary symplectic/anti-symplectic matrix 

v7 



F-iZ f) (198) 



It is readily verified that 
for all x, where 



f/lz^^f/l = , (199) 



7? + rr^ 

It is worth examining the map defined by Eqs. (|200p . (|20ip in a little more detail. 
In the first place observe that it is defined for all matrices F (including the singular 
ones), and not just for the matrices £ ESL(2,F(;). Furthermore, it is invertible, 
with inverse 

P = ^ '- h '- 203 

1 = \ — — (^04) 

77 + 77 

^ ^ 77(0 + 77-'&'^)+77-i(a'^ + 77^6) (2pg) 
77 + 77-1 

It is thus a bijection of the set of all 2 x 2 matrices over onto the set of all 
pairs G {¥^2)^. It obviously preserves the additive structure. Moreover, if we equip 
{¥cp)^ with the product rule 

(a, b) o (a', b') = {aa' + bb"^ , ab' + ba"^) (206) 

it also preserves the multiplicative structure (it is a ring isomorphism, in other 
words). Observe, further, that if F is any 2x2 matrix over F^, and if (a, 6) is the 
corresponding element of (F^) , then 

det F = a'^+i - b'^+^ (207) 

So we can identify ESL(2,Frf) with the set 

{(a, b) e (Frf2)2 : a-^+i - 6^+1 = ±1} (208) 

equipped with the above product rule. 

We are now ready to address the problem with which we started. For each 
{)<r<(f -2 define 

") (209) 
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(we define it like this, with a negative exponent of 9, with a view to later conve- 
nience). Wc then have 

FWr = Ur_logj(a+6?)-2'-) (210) 

So this gives us a natural way to represent ESL(2,F(i) as a group of permutation 
matrices. Note that the zero vector is excluded from the labelling scheme; however 
this does not matter because is invariant under the action of ESL(2,Frf). 
Finally observe that if 

-(;,-v-) 

then 

Aur = u^_i^ (212) 

So A is represented by a shift operator. This is a consequence of choosing the basis 
{77^^, 1} in Eq. (|193p . With a different choice of basis for F^2 one can arrange that 
any other order 2{d-\- 1) anti-symplectic becomes a shift operator. 

8. The MUB Cycling Problem: Preliminaries 

Wootters and Sussman [25] have shown that if d is a power of 2 there exists a 
Clifford unitary which cycles through a full set of mutually unbiassed bases (MUBs) . 
In the next two sections we address the problem of generalizing their result to the 
case when d is an odd prime power dimension. 

Let us stress that, throughout the remainder of this paper, we use such phrases 
as "the full set of MUBs" to mean "the full set of MUBs generated by acting on 
the standard basis with elements of the group f^" (which, as we will see, is just 
the set of MUBs orignally described by Wootters and Fields [33] )• It is possible to 
generate other MUBs by acting with the larger group £4. Also there is, of course, 
the possibility that there exist full sets of MUBs which are not covariant with 
respect to any version of the Clifford group. 

We will show: 

(1) There is no Clifford unitary which cycles through all the MUBs for any d. 

(2) If ci = 3 (mod 4) there is an anti-Clifford anti-unitary which cycles through 
all the MUBs. 

(3) For every d the MUBs can be divided into 2 groups of in such a way 
that there is a Clifford unitary which cycles through all the MUBs in each 
group separately. 

Let us stress once again that we are confining ourselves to the unitaries/anti- 
unitaries S f^. As Gross and Chaturvedi [5S] have observed the considerations 
which follow are not sufficient to exclude the possibility that there exist unitaries 
or (when d = 1 (mod 4)) anti- unitaries G £d which cycle through all the MUBs 
(except, of course, in prime dimension when the restricted and unrestricted groups 
are identical). 

In a subsequent paper this result will be used to show that minimum uncertainty 
states [251 I26[ 177] exist in every odd prime power dimension, thereby extending 
Wootters and Sussman's result [ISlllS], that such states exist in every even prime 
power dimension, and in a certain infinite class of dimensions = 3 (mod 4) . 

The purpose of this section is to fix notation, and to derive a formula which 
describes how an arbitrary matrix F G ESL(2,Fd) moves us between the different 
MUBs. In the next section we will prove our main result. 



29 



Let 



F 



a P 
7 S 



(213) 



be any element of SL(2,Fc;). It can be seen from Eq. (|76p that if /3 = then Up 
simply permutes and re-phases the standard basis. If, on the other hand, (3 =/= 
every element of Up has modulus . So if we define 




II ^ oo 



11 = oo 



(214) 



and 



(215) 



and allow fi to vary over the set F^; U {oo}, this will give us a full set of d + 1 
MUBs. The reader will easily perceive that, modulo some differences in the phase 
and labelling conventions, these are the MUBs originally described by Wootters 
and Fields [33 . 

Let us now consider the action of the extended Clifford group operators. For the 
displacement operators we have 



^tr((2a;-«2)«l)|^^2,_U2) 



/i 7^ 00 
fl — 00 



Turning to the case of the symplectic and anti-symplectic operators let 
be any element of ESL(2,F<j), and let A = detF. Then 



FH„ 



2 afj.+f3 \ / A 



.0 







' r 

^0 1^ 

.-1 0, 



1 / \7 TA^ + ^y 
' '-7 \ 



a 

( J 



(6 0' 



^ 7^ 00, 7/i + (5 7^ 
/i 7^ 00, 7/i + (5 = 
/i = 00, 7 7^ 
^ = 00, 7 = 



(216) 



(217) 



(218) 
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Consequently 



UF\^i,x) = < 



^ ^ (^) -"^^^ ^> ^ oo, + <5 ^ 

I (-7) f("T^') |oo, -7x) ^ 7^ 00, 7/1 + (5 = 



7 



7 ' 7 



(219) 

= 00, 7 7^ 



In particular the Wootters-Fields MUBs are the only MUBs which can be got by 
acting on the standard basis with the unitaries/anti-unitaries e 

9. The MUB Cycling Problem: Main Results 

We are now ready to prove our main results. The matrix is order p if /i 7^ 00, 
and order 4 if /i = 00. So if one wants to generate all the MUBs by acting on the 
standard basis with powers of the operators Uh^ one needs to use a minimum of 
pTi-i _|_ Y different operators. We will now show that it is possible to generate them 
using a much smaller number of operators. Specifically, we will show that in every 
odd prime power dimension it is possible to generate every MUB by acting on two 
distinguished bases with powers of a single unitary operator, and that if d = 3 
(mod 4) one can generate every MUB by acting on any single basis with powers of 
a single anti- unitary operator. 

It can be seen from Eq. (j216p that the displacement operators simply permute 
the vectors within a basis. So if one is interested in finding operators which cycle 
through the largest possible number of different bases one may confine oneself to 
the symplectic unitaries and anti-symplectic anti-unitaries. 

It follows from Eq. ((^1^)) that, for any F e ESL(2,Frf), 

UF\fi,x) - e'''-(^^--^fpUi),gpifi,x)) (220) 

for suitable functions Op, fp, gp- We will say that is a cycling matrix if one gets 
the complete set of MUBs by acting on any single MUB with powers of Up ■ 

Let us begin by observing that if is a cycling matrix, then so is every conjugate 
of F. In fact suppose that F is cycling, and suppose F' = SFS^^ for some S G 
ESL(2,Fd). The fact that F is cycling means that, for all /i, f™{y) = /i if and only 
if m = (mod d+ 1). Since f^i{n) = if and only if /]?(/s(m)) = /s(m) it follows 
that fpi has the same property, implying that F' is also a cycling matrix. 

It will be convenient to define the cycling index cp to the smallest positive integer 
m such that /™(0) = 0. In view of Eq. ([75]) cf is the smallest integer m such that 
F"^ is of the form 

2) (221) 

Since cp divides the order of F, and since the order of every type 1 matrix is 
< d — 1, we can immediately conclude that no type 1 matrix is a cycling matrix. 
It is also easily from Tables [1] and [2] together with Eq. (I158P that cf = 1 or p for 
every type 3 matrix. So if cycling matrices exist they have to be type 2. Every 
type 2 matrix is conjugate to a power of the anti-symplectic matrix 
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It is easily seen that if is a cycling matrix then so is A. We conclude that cycling 
matrices exist if and only if A is a cycling matrix. 
To calculate ca we use the fact 



^ = ^(n (223) 



where 



Consequently 



where 



S=\ (224) 



f3r = ^^^^Jr^ (226) 

(where we have used the fact that r/'^ = —77^^). So ca is the smallest positive 
integer m such that 

^m(d-l) ^ (227) 

Now /^"'(''^i) = 1 if and only if m{d — 1) = (mod 2{d + 1)). So ca is the smallest 
positive integer m such that 

m(d- 1) = 2fc(d+ 1) (228) 

for some k. 

Suppose that d = 4u + 1 for some u. Then Eq. (|228p becomes 

mu = k{2u + 1) (229) 

Since u and 2m + 1 are relatively prime we conclude that ca = ■ So A is not a 
cycling matrix. 

Suppose, on the other hand, that d = iu + 3 (mod 4). Then Eq. (|228p becomes 

m(2u+ 1) = 4fc(M+ 1) (230) 

Since 2u + 1 and 4(?i + 1) are relatively prime we conclude ca = d + 1. So A is a 
cycling matrix. 

Finally, note that if A is a cycling matrix then is also a cycling matrix if and 
only if r and d + 1 are relatively prime. In particular, the fact that d + 1 is even 
means that even powers of A are never cycling matrices. So there are no symplectic 
cycling matrices. 

We have thus proved 

Theorem 7. Let d be any odd prime power dimension. Then 

(1) There are no symplectic cycling matrices. 

(2) If d = 1 (mod 4) there are no anti- symplectic cycling matrices either. 

(3) If d = 3 (mod 4) an anti- symplectic matrix F is a cycling matrix if and 
only ifTT{F) ^ rf — ri~^ for some r relatively prime to d+ 1. 

At this point let us stress once again the point we made earlier, that we are only 
considering the operators in As Gross and Chaturvedi [86] have observed, the 
above considerations do not exclude the possibility that for some or all values of 
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d there is a cycling unitary G C(j; or that for some or all values of d = 1 (mod 4) 
there is a cycling anti- unitary g Ed- 

We now turn our attention to what, in the absence of a cycling matrix, might 
be considered the "next best thing" : namely, a matrix F with the property that 
the MUBs split into two groups of , in such a way that Up cycles through each 
group separately. We will refer to such matrices as half-cycling matrices. As with 
cycling matrices, if F is half-cycling, then so is every conjugate of F. 

For given F and /i define mi?(/i) to be the smallest positive integer m such that 
Ppip) = ^i- Clearly F is a half-cycling matrix if and only if mp{^) = ^^-i for all ^i. 

Since mp{^) is a factor of ord(_F) for all ^ no type 1 matrix can be a half-cycling 
matrix. Since mp{Q) — 1 or p for all type 3 matrices no type 3 matrix can be a 
half-cycling matrix. So, as with cycling matrices, if half-cycling matrices exist at 
all they must be type 2. It follows that every half-cycling matrix must be conjugate 
to a power of the matrix A defined in Eq. (|222p . 



To show that is a half-cycling matrix it is enough to show that rn^r (/i) 



d+l 
2 

for more than half the possible values of /i. So in the following we will assume that 
/i 7^ oo. Observe that if f^r{p-) = /i we must have 

iUAr''\fi,x) =e''\^i,kx) (231) 

for some 9, k. Equivalently 

UA^r\0,x) = e'^\0,kx) (232) 

where 

= H-^A^'^'H^ 

(/3mr— 1 f-^Pmr^ ( Pmrf-^ ' + iPmr-l - /Smr+l)/! + I3,nr)\ (^233) 



Referring to Eq. ((76)) we see that this implies 

PmrfJ-'^ — iPmr-l — /3mr+l)/^ " Pmr = (234) 

This is a quadratic equation in /i with discriminant 

iPmr-l - P,nr+lf + 4/3,^ = (jy™' - (- l)™'^,?"'" ^)2 (235) 

If the equation has solutions in the discriminant must e Q U {0}, implying 

which means that 77""' - (-1)™''^-™'- ^ if and only if 77"'^ - (_i)""-^-""- ^ 0. 
So f^r (/x) — fjL implies 

^2mr ^ (_;L)mr (237) 



or 



mr(d -1) = (mod2(d+l)) (238) 
Suppose, now, that d= \ (mod 4). Then Eq. (|238p reads 

mr I '^—r^ = ^mod ^^"^^ (239) 
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Since is relatively prime to this implies 

mr = ^mod (^40) 

We have thus shown that a necessary condition for f™r{fi) — ii is that mr is a 
multiple of But it is readily verified that 

k(d + l) 



fc(d+l) 



for all /i. So the condition is not only necessary but also sufficient. It follows that, 
for all 7^ oo, 

d+1 

"^^'■(/^)-|;ri±ij (242) 

In particular, is a half-cycling matrix if and only if r is relatively prime to 
Suppose, on the other hand, that d = 3 (mod 4). Then Eq. (|238p reads 

^ =0 (mod d+ 1) (243) 

So a necessary condition for /^-(/^) ~ fi is that mr is a multiple of rf + 1. But it is 
readily verified that 

K''^'' = (244) 

irrespective of the value of /i. So the condition is also sufficient. It follows that, for 
all ^ oo, 

"^^'-(^)-r^jrTT (245) 

[r, d+l\ 

In particular, is a half-cycling matrix if and only if r is an even integer such 
that I is relatively prime to 
We have thus proved 

Theorem 8. Let d be any odd prime power dimension. Then 

(a) A symplectic matrix F is a half-cycling matrix if and only if Tr(F) = 
^2r Jqj, some integer r relatively prime to 

(b) If d = 1 (mod A) an anti- symplectic matrix A is half-cycling if and only if 
Tt(F) = if — rj^^ for some odd integer r relatively prime to d + 1. 

(c) If d = 3 (mod A) there are no half-cycling antisymplectic matrices. 

10. Alternative Labelling Scheme 

If we label the MUBs in the way we have been doing so far, in terms of the 
parameter /i, the cycling structure described in the last section is not very obvious. 
We would like to construct a labelling scheme which makes it more explicit. 

This is, in essence, an easy problem to solve: all we have to do is choose an 
order 2{d + 1) anti-symplectic A and then label an arbitrary MUB by the number 
of applications of Ua required to reach it, starting from some distinguished MUB 
on the same orbit. If d = 3 (mod 4) it becomes easier still: all one has to do 
is count the number of applications of Ua needed to reach the MUB of interest 
starting from the standard basis. However, there is a slight complication due to 
the connection between the MUBs and the displacement operators: each MUB 
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determines a maximal commuting set of displacement operators and conversely 
(namely, the set of displacement operators which it diagonalizes). We would like 
our labelling scheme to make this relationship transparent. To do that we will make 
use of the labelling of the displacement operators described in Section [T] 

It is easily seen that the maximal commuting sets of displacement operators 
corresponding to the Wootters-Fields MUBs are all of the form 

{D^u-weFd} (246) 

for some fixed u 7^ 0. In terms of the labelling of the displacement operators 
described in Section [7] we have 

WUr ^ Ur-(d+l)\oggW (247) 

for all non-zero w G F^. This suggests replacing the single integer r with a pair 
of indices s, t, defined so that s (respectively t) is the unique integer in the range 
< s < d (respectively < t < d — 2) such that 

r^s + t{d+l) (248) 

(in other words t and s are, respectively, the quotient and remainder of r on division 
by d + 1). In terms of this double index notation we have 

wus^t ^ Usj-ioggw (249) 

for all non-zero w €¥d. So the maximal sets of commuting displacement operators 
corresponding to the Wootters-Fields MUBs are the d+ 1 sets 

Ss = {Du,t:0<t<d-2}U{Do} (250) 

obtained as s ranges over the interval < s < d. Thus, the index s labels the 
maximal set and the index t labels the individual displacement operators within 
the set. 

We now define, corresponding to each Sg, the d operators 

l + ^'c.-*'(''"*-)i?„^ J (251) 
t=o / 

It is readily confirmed that the Ps ^ are rank 1 projection operators and, further- 
more, that 

tr(P,,,P,^,0 = i 'Z''^ (252) 

So if we set 

Ps,x = \s,x)c c{s,x\ (253) 
("c" for "cycling") the vectors \s,x)c will constitute the full set of MUBs, with s 
labelling the basis and x labelling the individual vector within the basis. Moreover, 
{\s, x) : X G ¥d} is the basis which diagonalizes the displacement operators in Ss- 

To find the relation between this labelling scheme and the one used previously 
observe that, in view of Eq. ()216|) . 

= "^^^ (254) 
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implying 

\s,x)cMtJ-s,Kx) (255) 
where the notation "=" means "equals up to a phase" , and where 

v'^'+v-'-' ^ 2 (256) 

X, = \ vS-^+v-'S-^^ ^ (257) 



^1 

We also have the inverse relation 



2 



A 

where 



\^^:X)^S^,,^x)c (258) 



d- l 
. 1 



(259) 



H = oo 



In particular |0,a;)c is the standard basis. 

To calculate the action of Up, for arbitrary F e ESL{2,¥d), let /f(s) be the 
unique integer in the range < /i?(s) < d and gris) the unique non-zero element 
of ¥d such that 

s ~ loggiap - bprr^n = /f(s) + (d + 1) log, (^^(s)) (260) 

(so logg {gpis)) is the quotient and /f(s) is the remainder on division by d + 1). 
Then 

^Ps,t = U/^(s),t+loge(gF(s)) (261) 

In view of Eq. (|25ip this means 

C/fP... - ^ (E^""^^'"*^''^^^^^^u,,u,.) (262) 

implying 

f/F|s,a;)c=|/F(s),ffF(s)x)c (263) 
Finally, let us consider the action of the operators Ua, Ua2. We have a a = rj, 
Ba = 0. In view of Eq. (|206p this means a^m = rj™, = for all m. So 

Ua2\s,x)c=1 ' (264) 

I |s + 2, a::)c s = d — 1 or d 

(where addition in the first argument is mod d + 1). We see that, for all d, Ua^ 
cycles through the even-index and odd-index MUBs separately. 
Turning to Ua we have 

If d = 1 (mod 4) then is even, so Ua, hke [/^2, cycles through the even- index 
and odd-index MUBs separately. But if d = 3 (mod 4) then is an odd integer 
relatively prime to d+ 1, so Ua cycles through all the MUBs. 
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Lastly, let us note that, although A is a eycling matrix when d — 3 (mod 4), it 
is rather more convenient to use as this increases the index by 1 each time: 

' < ' (266) 
s = a 

11. Concluding Remark 

In this paper we have confined ourselves to the group It would obviously 
be interesting to see how far our results extend to the larger group £d- It would 
also be interesting to examine the case of even prime power dimension. We hope 
to address these questions in future publications. 
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